What we (don't) know about black hole formation in high-energy collisions 
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Higher-dimensional scenarios allow for the formation of mini-black holes from TeV-scale particle 
collisions. The purpose of this paper is to review and compare different methods for the estimate 
of the total gravitational energy emitted in this process. To date, black hole formation has mainly 
been studied using an apparent horizon search technique. This approach yields only an upper 
bound on the gravitational energy emitted during black hole formation. Alternative calculations 
based on instantaneous collisions of point particles and black hole perturbation theory suggest that 
the emitted gravitational energy may be smaller. New and more refined methods may be necessary 
to accurately describe black hole formation in high-energy particle collisions. 
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I. INTRODUCTION 

The standard model of particle physics has been suc- 
cessfully tested up to energies of ~ 1 TeV. However, 
its foundations are still mysterious. Much effort has 
been spent to explain the hierarchy problem, i.e. the 
huge difference between the electroweak scale, ttiew ~ 
300 GeV, and the Planck scale, M P1 - 10 19 GeV. While 
electroweak interactions have been probed at distances 
m EW ~ 10~ 16 cm, gravitational forces have not been 
probed at distances Mp^ ~ 10~ 33 cm. Gravity has only 
been accurately measured in the ~ 0.01 cm range [1]. 

If gravity is modified at scales smaller than 1 mm, the 
hierarchy problem can be solved by assuming the exis- 
tence of n compact extra dimensions of length ~ R [2]. 
Gauss's law in D = 4+n dimensions implies that two test 
masses mi, mi at a distance r <^ R feel a gravitational 
potential 

V (r) = G D ^^, r«R, (1) 

where Go is the D-dimensional Newton constant. The 
four-dimensional Newtonian potential is recovered at dis- 
tances r » R, where r" ~ R n and G 4 = G D /R n . The 
D-dimensional Planck mass, Mp^n, is obtained by equat- 
ing the Schwarzschild radius of an object of mass m to 
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its Compton wavelength, A = 1/m: 

(%d) D - 2 ~^. (2) 

(Here and throughout the paper we set h = 1 and c = 1.) 
Mpi o is related to the effective four-dimensional Planck 
scale Mpi by M|j ~ {Mv\,T>) D ^ 2 R n . The hierarchy prob- 
lem is solved by imposing equal scales in the higher- 
dimensional setting, i.e., ttiew — Mpi,D< This condition 
relates the size R of the extra dimensions to the number 
of extra dimensions n. A single extra dimension implies 
deviations from Newtonian gravity over solar system dis- 
tances. It is thus excluded empirically. If n = 2, the size 
of the extra dimensions is R ~ 0.3 mm. This value has 
been recently ruled out by experiments with torsion pen- 
dulums [1]. However, any n > 3 gives modifications of 
Newtonian gravity at distances smaller than those cur- 
rently probed by experiment. 

An exciting consequence of TeV-scale gravity is the 
possibility of production of black holes (BHs) in particle 
colliders [3, 4] and ultra high energy cosmic ray interac- 
tions with the atmosphere [5]. A naive estimate of the 
cross section for Mp^d ~ 1 TeV predicts that super- TeV 
particle colliders will produce BHs at a rate of few per 
second. 

The BH starts to decay after forming. First, it radiates 
all the excess multipole moments. The BH has initially 
a highly asymmetric shape. It then settles down to a 
stationary state. (In four dimensions it must be a Kerr- 
Newman BH, by the uniqueness theorem.) This phase 
is commonly called balding phase. The endpoint of the 
balding phase is a spinning BH. This BH starts to lose an- 
gular momentum (spin-down phase). Page [6] has shown 
that a spinning BH radiates mainly along the equatorial 



2 



plane, and that this radiation carries away most of the 
BH angular momentum. After completing the spin-down 
phase, the BH radiates its remaining degrees of freedom 
through Hawking radiation (Schwarzschild phase) [7, 8]. 
The endpoint of the Schwarzschild phase is either explo- 
sive or a BH remnant is left [9]. 

The focus of this paper is on the formation of a BH 
in high-energy particle collisions. In particular, we are 
interested in the total gravitational energy radiated. Dif- 
ferent methods have been discussed in the literature. The 
standard procedure is to describe the incoming particles 
by two Aichclburg-Sexl shock waves [10], and then find a 
closed trapped surface in the union of these shock waves 
[11-18]. The shock wave approach can be used to find 
lower bounds on the mass of the newly formed BH, i.e., 
upper bounds on the amount of gravitational radiation 
emitted during the process. However, there is little dis- 
cussion in the literature on how well these results approx- 
imate the actual gravitational emission. The aim of this 
paper is to provide a critical discussion of the different 
methods. In Section II we describe briefly the standard 
Aichelburg-Sexl shock wave technique, the main results 
of this approach, and its shortcomings. In Section III 
we discuss two alternative techniques. Conclusions are 
presented in Section IV. 



II. BLACK HOLE FORMATION VIA 
AICHELBURG-SEXL SHOCK WAVES 

Consider a D-dimensional Schwarzschild- Tanghcrlini 
solution with mass M [19]. The metric in spherical co- 
ordinates is 



ds 2 = - 1 - 



IQttGdM 



1 



(D - 2)f! D _ 2 r D ~ 3 
-l 



dt 2 



1 - 



16ttG d M 



1 



{D-2)n D _ 2 r D ~ 3 



dr 2 + r 2 dn 2 D _ 2 (3) 



where is the volume of the unit (D — 2)-dimensional 

sphere 



D-2 



2tt( d - 1 )/ 2 

r[(i>-i)/2]' 



(4) 



and T[z] is Euler's Gamma function. Boosting this solu- 
tion to very large values of the Lorentz factor 7 at fixed 
total energy fj, = M7, we obtain the Aichclburg-Sexl so- 
lution 

ds 2 = -dudv + (dx 1 ) 2 + <5>{x l )8{u)du 2 , (5) 
where $ depends only on the transverse radius p 2 — x l Xi\ 



$ = -8G D n\ogp , D = 4 
WttGdP 



$ = 



Q D -s(D - 4)pC-4 



, D >4. 



(6) 



Equation (5) describes a particle with energy fi moving in 
the +z-direction. The spacetime is flat outside the null 



TABLE I: Upper limits on the efficiency of gravitational ra- 
diation emission for different spacetime dimensions D, using 
the trapped surface method. 
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plane u = 0. If we consider an identical shock wave trav- 
eling along v = in the — z direction, the two solutions 
can be superposed to yield a solution for two colliding 
shock waves. BH formation can be studied by identify- 
ing a future trapped surface in the solution, with no need 
to calculate the gravitational field [11-18, 20]. If the col- 
lision is head-on, it is easy to find a trapped surface in 
the union of two flat disks with radii 



\ ^D-3 



l/(D-3) 



(7) 



Since the BH horizon is always in the exterior region, the 
method gives a lower bound on the final BH mass. If the 
final BH is non-spinning, its mass is 



Mbh > M 



(D-2)n D - 2 



2n 



D-3 



(8) 



The total gravitational energy radiated is 2/i — Mbh , an d 
the efficiency is e = 1 — Mbh/^H- An efficiency of 0% 
implies that no gravitational radiation is emitted. An 
upper limit on the efficiency follows from Eq. (8). Some 
values are shown in Table I. 

The above formalism can be extended to the study 
of collisions with nonzero impact parameter [13-15, 17]. 
This makes possible the computation of the BH produc- 
tion cross-section. 



A. Accuracy of the results 

The trapped surface method gives only a lower bound 
on the BH mass. The BH mass can be anything between 
this bound and the center of mass (cm.) energy of the 
collision. D'Eath and Payne studied this problem in more 
detail for the four-dimensional case [12]. Bondi's news 
function, which describes the emission of gravitational 
radiation, can be written as an infinite series around the 
collision axis. The first term of the series is Eq. (8). 
Making some assumptions on the angular dependence of 
the radiation, and extrapolating off the axis, the second 
term is found to decrease the efficiency for gravitational 
wave generation in head-on collisions to 16%. Although 
this derivation relies on various approximations, the re- 
duction in efficiency is significant, and may signal that 
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the series is slowly converging. The situation seems to 
get even worse in higher dimensions (see Section III) and 
for large impact parameters. As an illustration of how 
the trapped surface method may lead to an inaccurate 
estimate of the BH mass, let us consider the collision 
of two non-spinning BHs initially at rest. Using simi- 
lar arguments to the trapped surface method, Hawking 
[21] placed an upper limit of 29.3% on the efficiency of 
gravitational wave generation. The exact result can be 
obtained through a numerical solution of Einstein's equa- 
tions [22-28], and is around 0.1%, i.e. two orders of mag- 
nitude smaller than the Hawking bound. The results for 
ultra-relativistic collisions are likely to be more accurate. 
However, the total energy radiated may still be much 
smaller than the upper limits of Table I. 



B. Finite size particles 

The Aichelburg-Sexl solution describes the metric of a 
massless pointlike particle with a very large boost. Clas- 
sically, the point-particle assumption is accurate along 
the collision axis because of the large Lorentz contrac- 
tion due to the boost. However, it fails for directions 
transversal to the motion. If the colliding particles are 
strings, string size effects can be modeled by consider- 
ing beam-beam collisions, with beam sizes of order As, 
where As > Ipi [20]. The effects of finite-size transver- 
sal dimensions have been studied analytically by Kohl- 
prath and Veneziano [20], and lead to a smaller cross 
section. (We are not aware of any numerical study.) For 
a reasonable cross section the cm. energy should satisfy 
E > M Pl {\ s /l Pl ) D -\ 



C. Spin effects 

In general, colliding particles have intrinsic spin and 
should be modeled by metrics other than the Tanghcr- 
lini metric. A naive generalization of the Aichelburg-Sexl 
approach is to boost the rotating Myers-Perry metric 
[29]. This has been done by several authors [30]. The 
results are cumbersome enough to make the trapped sur- 
face method very difficult to implement. A seemingly 
better candidate to model spinning high-energy particles 
was recently proposed in Ref. [31]. This model consists in 
a spinning radiation beam-pulse which includes dragging 
effects, in contrast to the boosted Kerr metric. 

Another question concerns the angular momentum of 
the BH, given the intrinsic spin and angular momentum 
of the colliding particles. (This is also an important point 
in astrophysics [32].) The cross section and the angular 
momentum of the final BH can be estimated by assum- 
ing the net angular momentum carried by gravitational 
waves to be negligible [8, 13]. However, a more accurate 
investigation of angular momentum effects is needed. 



D. Charge effects 

String length and spin effects should be important if 
the incoming particles have energy close to the Planck 
scale. It is likely that both of these effects are suppressed 
for super-Planckian energies. Charge effects are expected 
to dominate at very high energies because gauge fields are 
confined on the brane and decay more slowly than the 
gravitational field of a neutral particle [33] . Estimates of 
charge effects in BH production from high-energy colli- 
sions are not yet available in the literature. A first at- 
tempt was presented in Ref. [34] using a perturbative 
method. 



III. OTHER METHODS TO ESTIMATE THE 
ENERGY RELEASED 

Because of the problems listed above, it is desirable to 
explore different methods. In this section we discuss two 
possible approaches. Although they are not free from 
shortcomings, quantitative results on the BH formation 
process can be obtained. These techniques give consis- 
tent results in D = 4. 

The formalism to handle gravitational waves in D- 
dimensional flat spacetimes is discussed in Ref. [35] . The 
nonlincarity of Einstein's equations makes the treatment 
of the gravitational radiation problem difficult. It is a 
standard procedure to work only with the weak radia- 
tive solution. In this limit, the energy-momentum self- 
interaction term of the gravitational wave can be ne- 
glected. This approach is justified in most problems 
where the total amount of gravitational radiation re- 
leased is negligible in comparison to the total energy con- 
tent of the spacetime. (High-energy collisions of two BHs 
require in principle the inclusion of nonlinear effects.) Let 
us assume an asymptotically flat D-dimensional space- 
time with metric <? M „ = r\^ v + h^ u . At this linearized 
level it can be shown that 

(i) gravitational waves in a D-dimensional spacetime 
have D(D — 3)/2 independent polarizations; 

(ii) the perturbation amplitude in the wave zone, h^, 
can be expressed in terms of the energy momentum ten- 
sor as 



I 



( 27rr )(U-2)/2 « 

tp-VS^t-lx-x'l.xO 



(9) 



/ D-4 \ 

where d + 2 stands for the ^^th derivative with re- 



spect to time and 



D-2 



(10) 



The Fourier transform and the energy spectrum are (see 
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Ref. [35] for details) 

(11) 

(T^( w ,k)r;(,,k) - _l_|T\( w ,k)| 2 ) ,(i2) 

respectively. 



TABLE II: Efficiency for gravitational radiation generation in 
head-on collisions for different spacetime dimensions D, using 
the instantaneous collision method. 
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A. Instantaneous Collisions in Even .D-Dimensions 

In general, two scattering bodies release gravitational 
energy due to momentum exchange. If the collision is 
hard, i.e. the incoming and outgoing trajectories have 
constant velocities, the metric perturbation and the re- 
leased energy can be computed exactly. (This method 
was first derived by Weinberg [36, 37] and later explored 
in Ref. [38].) These calculations assume an instantaneous 
collision, and are valid for arbitrary velocities and low en- 
ergies. The resulting spectrum is flat in four dimensions 
[39], and a cutoff frequency is needed to obtain a finite 
total energy. A suitable cutoff enables us to estimate 
the total energy radiated by the collision of two ultra- 
relativistic BHs. This method has been recently general- 
ized to higher dimensions in Ref. [35] . In what follows we 
revisit this result and estimate the total energy radiated 
in the high-energy collision of two D-dimensional BHs. 

Consider a system of freely moving particles with D- 
momenta P?, energies Ei and (D— l)-vclocities v. These 
quantities change abruptly at t = to corresponding 
primed quantities due to the collision. The energy- 
momentum tensor is 



Eq. (13) and (12), we find 



Ei 



_L_^-V-v't)e(i) 



(13) 



Substituting Eq. (13) in Eqs. (11) and (12), the metric 
perturbation h^ v and the radiation spectrum can be ob- 
tained. Let us consider a head-on collision of a particle 
with mass toi and Lorentz factor 71 with a particle of 
mass mi and Lorentz factor 72 in the cm. frame. We 
assume without loss of generality that the motion is in 
the (xd-i,%d) plane: 

Pi = 7imi(l, 0, 0, vi sin#i, v\ cos#i) , 

P[ = (^,0,0, ...,0,0), (14) 

P 2 = 72m2(l,0,0, -W2sin#i, -w 2 cos#i) , 

P' 2 = (£ 2 , 0,0, ...,0,0). (15) 

Momentum conservation leads to the additional relation 
^\m\V\ = 72m 2 t>2. Substituting Eqs. (14) and (15) in 



d 2 E 



2G D D -3"f1m 2 1 vl(v 1 +v 2 ) 2 sin 1 9 1 uj 



D-4 



dujdfl {2ir) D - 2 D-2 (l-«iCOS(9i) 2 (l + v 2 cos(9i) 2 ' 

(16) 

For arbitrary (even) dimensions D > 4 the spectrum is 
not flat. The integration of Eq. (16) from u = to a 
cutoff frequency uj c gives 



dE 



2G, 



1 7 2 to 2 tj 2 (tji + v 2 ) 2 sin 4 9\ u>, 



D-3 



dn (2tt) d - 2 D- 2 (1-viCOs6>i) 2 (1 + u 2 cos(9i) 2 ' 

(17) 

If a BH forms in the collision, the effective timescale for 
the process is r ~ r+. This suggests to approximate the 
cutoff frequency lu c by the inverse of the BH radius. For 
a collision between equal-mass particles (mi = m 2 = to, 
v\ = V2 = v) this yields a total energy 



E : 



2 3 - g 7 2 TO 2 
AiT 2 [(D- l)/2] 



where M ~ 27777. The efficiency E/(2my) is 



2 1 



-D 



T 2 [(D -l)/2] 



(18) 



(19) 



Values of e are listed in Table II for different dimensions. 
The results are in agreement with the four-dimensional 
estimate of D'Eath and Payne [12]. However, the total 
energy decreases with the spacetime dimension, in dis- 
agreement with the estimate of Ref. [13]. 

The instantaneous collision approach seems to suggest 
that the trapped surface method overestimates the total 
energy emitted in the collision. The approach described 
in this section can be generalized to describe the collision 
of rotating bodies. Extension to rotating BHs would also 
be of interest for the computation of gravitational radia- 
tion from gamma-ray bursts [40] . 



B. Perturbation around the black hole background 

The collision of a BH with an ultra-relativistic particle 
was studied in detail by Cardoso and Lemos [41] in four 
dimensions, and by Bcrti et al. for D > 4 [42]. In the per- 
turbative approach, the BH-particle system is described 
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TABLE III: Efficiency for gravitational radiation generation 
in head-on collisions for different spacetime dimensions D, 
using the perturbed Schwarzschild-Tangherlini BH method. 
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by a perturbed Schwarzschild-Tangherlini BH, where the 
perturbation h^ v is induced by the infalling particle. Ex- 
panding Einstein's equations to first order in h^ Vl the 
problem can be expressed as a second order differential 
equation for h^ v . Although the formalism is strictly valid 
only for a colliding particle with small energy E com- 
pared to the BH mass, the results can be extrapolated 
to E <~ Mbh [42] . The results of this analysis are given 
in Table III. They are in qualitative agreement with the 
predictions of the instantaneous collision approach. 

The total energy radiated by a freely-falling particle in 
the four-dimensional Schwarzschild spacetime [43] is in 
very good agreement with numerical simulations of BH 
head-on collisions. An extensive comparison of numeri- 
cal and perturbative results can be found in Refs. [22- 
25]. According to perturbation theory, the total en- 
ergy radiated by a test particle of mass m falling ra- 
dially from rest at infinity into a BH of mass M 3> m 
is E = 0.0104to 2 /M [43]. If m is replaced by the re- 
duced mass fi, this result agrees with the numerical re- 
sult E ~ 0.0013 (see Fig. 14 and Sec. IV of Rcf. [25]). 
The slight discrepancy between perturbative and numer- 
ical results can be quantitatively explained by consider- 
ing three fudge factors in the perturbative analysis [25]: 
F ra , which accounts for the finite initial infall distance in 
the numerical simulations; Fh, which accounts for tidal 
deformations heating up the BH horizon; .Fabs, which ac- 
counts for the reabsorption of the gravitational waves by 
the BH. In our case the process starts at infinite separa- 
tion, thus F ro = 1, F h ~ 0.86, and F ahs ^ 0.99. 

Perturbation theory in the close-limit approximation 
and numerical simulations are also consistent for four- 
dimensional ultra-relativistic collisions. Numerical stud- 
ies of boosted BHs with fixed initial separation show that 
the energy emission saturates at E ~ 0.0 1M for very 
large initial BH momenta (see, e.g., Fig. 2 of Rcf. [27]). 



Ref. [26] combines Newtonian dynamics and numerical 
simulations of boosted BHs to conclude that the maxi- 
mum energy emission could actually be much lower than 
the above value, E < 0.0016M. Four-dimensional inves- 
tigations also suggest that most of the radiation is emit- 
ted in the ringdown phase. This is confirmed by compar- 
ing perturbative calculations to post-Newtonian calcula- 
tions [44]. The bremsstrahlung radiation of a particle at 
distance larger than AM (in Schwarzschild coordinates) 
contributes only <~ 3 % of the total energy emitted. 

In conclusion, numerical results in four dimensions in- 
dicate that perturbation theory likely overestimates the 
emitted radiation. The generalization of these results to 
higher dimensions is nontrivial, and the values in Tables 
II and III could underestimate the efficiency for gravi- 
tational wave generation in high-energy collision of two 
equal mass particles. 



IV. CONCLUSIONS 

We reviewed and compared different approaches to 
computing the gravitational energy released during back 
hole formation in high-energy particle collisions. While 
in four dimensions there is good agreement between all 
these methods, results differ significantly in higher di- 
mensions. The straightforward conclusion is that new 
techniques and refinements in previous calculations are 
needed to obtain a quantitative understanding of BH for- 
mation in higher-dimensional spacetimes. This is partic- 
ularly important for high D, where different approaches 
do not yield consistent results, and for non head-on col- 
lisions, where the trapped surface method is expected to 
be less reliable. 
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